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The recent surge of interest in phase change materials GeTe, Ge 2 Sb 2 Tes, and related compounds 
motivated us to revisit the structural phase transition in GeTe in more details than was done before. 
Rhombohedral-to-cubic ferroelectric phase transition in GeTe has been studied by high resolution 
neutron powder diffraction on a spallation neutron source. We determined the temperature de¬ 
pendence of the structural parameters in a wide temperature range extending from 309 to 973 K. 
Results of our studies clearly show an anomalous volume contraction of 0.6% at the phase transition 
from the rhombohedral to cubic phase. In order to better understand the phase transition and the 
associated anomalous volume decrease in GeTe we have performed phonon calculations based on 
the density functional theory. Results of the present investigations are also discussed with respect 
to the experimental data obtained for single crystals of GeTe. 

PACS numbers: 


I. INTRODUCTION 

Phase change materials possess unique properties that 
hold considerable promise for applications in data stor¬ 
age ID- They can be rapidly and reversibly switched be¬ 
tween amorphous and crystalline states, which differ sub¬ 
stantially in their properties. Recently materials such as 
Ge 2 Sb 2 Te 5 and Ag and In doped Sb 2 Te have been discov¬ 
ered to crystallize rapidly enough to enable competitive 
solutions for rewritable optical data storage. 

The related binary material GeTe has drawn consid¬ 
erable interest because of its higher crystalline temper¬ 
atures and better data retention at high temperatures 
compared to Ge 2 Sb 2 Te 5 . Apart from the application in 
data storage, GeTe may have potential use as a ther¬ 
moelectric material [2- These have motivated us to re¬ 
visit the high temperature ferroelectric phase transition 
in GeTe. 

GeTe is a narrow band-gap semiconductor [3] and is 
ferroelectric at room temperature with a Curie temper¬ 
ature of about 705 K. The low temperature ferroelectric 
phase has a rhomboliedraly distorted NaCl type crystal 
structure with the space group i?3m I3HS1- Structural 
distortions involve relative displacement of the Ge and 
Te sublattices along the body cell diagonal and subse¬ 
quent rhombohedral shear deformation along [111] direc¬ 
tion which changes the rhombohedral angle from its fee 
value of 60° to a. The Ge and Te atoms are six fold 
coordinated by each other with three shorter (2.83 A) 
and three longer (3.15 A) bonds. This is often described 
as Peierls distortion due to reduced coupling between 
the p-type orbitals that constitute basis of the bonding 
in GeTe. 

GeTe undergoes a ferroelectric phase transition in 
which the low temperature rhombohedral R2>m structure 
transform to the cubic Fm3m structure at high temper¬ 
ature of about 600-700 K [7, j3[. The transition tem¬ 


perature T c depends on the sample stoichiometry and 
carrier concentration [5j. Ferroelectric phase transition 
in GeTe was considered to be displacive in its origin 
o mm mi- Recently, however, the displacive character 
of the rhombohedral-to-cubic phase transition in GeTe 
has been contested by Fons at al. m and Matsunaga 
et al. m- According to their studies the displacive na¬ 
ture of this phase transition was due to the misinterpre¬ 
tation of the Bragg diffraction results as the structure 
determination based only on the Bragg intensities gives 
information about the average structure, but not about 
the system local structure. Information about the lo¬ 
cal structure can be obtained from the total scattering 
data including the diffuse scattering up to a very high 
Q value and from the pair-distribution function (PDF) 
analysis (14] . Such investigations has been performed 
on the X-ray diffraction data by Matsunaga et al. [13] 
and lead to the conclusion that in the local scale the 
high temperature phase of GeTe still exhibits distinct 
short and long Ge-Te bonds, contrary to the conven¬ 
tional structure refinement with only the Bragg intensi¬ 
ties considered which suggests that the high temperature 
cubic phase reveals solely one type of the Ge-Te bond. 
Additionally, the two distinct bond distances observed in 
the local scale hardly change across the rhombohedral-to- 
cubic phase transition. Results of Matsunaga et al. E3 
were also supported by the EXAFS studies of Fons et al. 
m- These findings suggested that the phase transition 
in GeTe is not displacive but order-disorder type. On 
the other hand, the most recent lattice dynamical cal¬ 
culations m based on density functional theory (DFT) 
show that the rhombohedral-to-cubic phase transition in 
GeTe is indeed displacive in its origin and becomes driven 
by the condensation of exactly three components of the 
triply degenerate optical transverse soft phonon mode at 
the Brillouin zone center. Moreover, the displacive char¬ 
acter of the phase transition in GeTe has been further 
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supported by the recent electron and X-ray diffraction 
studies as well as the Raman scattering experiments of 
Polking et al. m- 

We have revisited the ferroelectric phase transition in 
GeTe by the high resolution neutron powder diffraction 
on a modern high power spallation neutron source and 
determined the temperature variation of the lattice pa¬ 
rameters, unit cell volume, positional parameters, and 
bond distances more accurately and in much finer tem¬ 
perature steps across the rhombohedral-to-cubic phase 
transition than were done for this compound before [3]. 
Also, the DFT phonon calculations were performed to 
better understand the phase transition and the associ¬ 
ated anomalous volume contraction in GeTe. 


II. EXPERIMENTAL AND CALCULATION 
METHODS 

GeTe powder samples were obtained from Alfa Aesar. 
The samples are claimed to be 99.999% pure and are 
200 Mesh. We checked the samples by the X-ray pow¬ 
der diffraction and found the presence of small amount 
of GeOo impurity. Neutron powder diffraction (NPD) 
measurements were performed on the time-of-flight pow¬ 
der diffractometer, POWGEN, located at the Spallation 
Neutron Source at Oak Ridge National Laboratory. The 
data were collected with neutrons of central wavelength 
1.333 A, covering d-spacing range from 0.42 to 5.4 A. 
Approximately 5 g of GeTe sample was loaded in a vana¬ 
dium container of 10 mm diameter and measured in a 
traditional ILL furnace within the temperature range of 
309-973 K. Structure refinement was carried out using 
the FullProf suite [16l . 

The present theoretical studies use the DFT method 
implemented in the VASP code m together with the 
direct method approach [T8] to provide necessary in¬ 
put data to calculate the temperature dependence of the 
mean-squared vibrational amplitudes ( U-ij) of the Ge and 
Te atoms in the low and high temperature phases of 
GeTe within the harmonic theory. The Uij tensor is ob¬ 
tained from the calculated diagonal and off-diagonal par¬ 
tial phonon densities of states |H] ■ The volume thermal 
expansions for the low and high temperature structures 
of GeTe are evaluated according the quasiharmonic ap¬ 
proximation (QHA) [20]. Details of the lattice dynamics 
calculations on GeTe system can be found in Ref. m- 
The pair-distribution functions of GeTe were obtained 
using the PDFgui program [2T :. 

III. EXPERIMENTAL RESULTS 
A. Crystal structure and unit cell transformation 

Before we discuss the present NPD results and com¬ 
pare them with the results of neutron diffraction exper¬ 
iments on single crystals [3], we describe the different 
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FIG. 1: (Color online) Structure of the low temperature 
phase of GeTe (space group R3m) shown in the pseudo-cubic 
(thick black lines), hexagonal (dotted lines), and rhombohe- 
dral (thick blue/gray lines) representations. 


units cells used in the two investigations as well as the 
relationship between them. The low temperature struc¬ 
ture of GeTe with the space group R2>m can be expressed 
in the pseudo-cubic, rhombohedral or hexagonal crystal¬ 
lographic representations as illustrated in Fig. [I] The 
distorted rocksalt structure with the lattice parameter a c 
and angle a is related to hexagonal unit cell with the lat¬ 
tice constants a and c via the relation a = 2a c sin(a/2) 
and c = a c y/3 + 6 cos a [22]. In the hexagonal repre¬ 
sentation the Ge and Te atoms occupy 3a(0,0,:r) and 
3a(0,0, 1 — x) Wyckoff positions, respectively. In addi¬ 
tion, the distortion parameter Aa: = 0.25 — XQ e = ^Te — 
0.75 describes the relative shift of the Ge and Te sublat¬ 
tices from the values xc e = 0.25 and xtc = 0.75 which 
are characteristic for the cubic GeTe structure. The de¬ 
gree of distortion from the cubic NaCl-type structure is 
also reflected in the deviation Aa of the rhombohedral 
angle a from the cubic value of 90° (Aa = 90° — a). 

We recall that Chattopadhyay et al. [8] performed 
structural refinement of their single crystal data using 
the pseudo-cubic cell with the F13m setting whereas we 
have done structure refinement of our NPD data in the 
hexagonal setting of a rhombohedral cell. At 309 K, 
the hexagonal cell is described by the lattice parame¬ 
ters ah = 4.1651 A and Ch = 10.6704 A, correspond¬ 
ing to the pseudo-cubic cell with the lattice parameters 
a c = 5.9818 A and a c = 88.2615°. Although we re¬ 
fined the powder diffraction data in the hexagonal cell 
we choose the pseudo-cubic representation as the most 
convenient to make comparison between the present re¬ 
sults and those of Chattopadhyay et al. [8]. The usage 
of the pseudo-cubic representations is also advantageous 
because a = 90° in the high temperature phase. 
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FIG. 2: (Color online) Temperature dependence of diffraction 
diagrams for GeTe measured at wide d-range (a) and d-ranges 
corresponding respectively to the cubic (222) c (b), (220) c (c), 
and (200) c (d) reflections. Note the phase transition at T c = 
600 K, development of the rhombohedral distortions below T c , 
and vanishing splitting of the diffraction peaks at the onset 
of the rhombohedral-to-cubic phase transition. 
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FIG. 3: (Color online) Results of the Rietveld profile refine¬ 
ment for the low temperature rhombohedral phase of GeTe 
(a) and the high temperature cubic phase of GeTe (b). The 
refinements were performed by taking into account the impu¬ 
rity phase GeC> 2 . The orange and green small vertical lines 
indicate the positions of the Bragg peaks corresponding to 
the GeTe and GeC >2 phases, respectively. The blue curves 
at the bottom denote differences between the measured and 
calculated intensities. The conventional Rietveld discrepancy 
parameters for the refinement (a) were R p = 20.9%, R wp = 
21.3%, Re = 8.48% and y 2 = 6.323 and those for the refine¬ 
ment (b) were R p = 26.6%, R wp = 18.9%, R e = 11.9% and 
X 2 = 2.509. 


B. Neutron Powder Diffraction Results 

Figures[2](a)-(d) show the temperature variation of the 
diffraction diagram represented by contour plots. The 
large d-range is presented in Fig.[2|a), while the d-ranges 
corresponding to the cubic (222) c , (220) c , and (200) c 
reflections are depicted in Fig. [2j)c)-(d). These reflec¬ 
tions are split due to distortions of the Peierls type be¬ 
low T c = 600 K. The splitting progressively diminishes 
with increasing temperature and finally the double-peak 
structure disappears at the onset of structural transfor¬ 
mation. The single peak observed above 600 K is a direct 
evidence of the cubic symmetry of GeTe. 

Diffraction intensities of the GeTe samples have been 
refined together with the impurity Ge02 phase. Re¬ 
sults of the Rietveld refinements performed at 310 K 
(low temperature structure) and 923 K (high tempera¬ 
ture structure) are shown in Fig. [3j The conventional Ri¬ 
etveld discrepancy parameters for the refinement (a) were 
R p = 20.9 %,R wp = 21.3 %,R e = 8.48% and x 2 = 6.323 
and those for the refinement (b) were R p = 26.6%, R wp = 
18.9%, R e = 11.9% and y 2 = 2.509. The suffixes p, wp 
and e mean profile, weighted profile and expected from 
the counting statistics, respectively. We preferred to call 


these parameters discrepancy rather than agreement pa¬ 
rameters because these parameters are larger when the 
discrepancy (and not the agreement) is larger. It occurs 
that the impurity phase does not affect the refined pa¬ 
rameters neither of the low nor high temperature GeTe 
phases. We performed such refinements for the data mea¬ 
sured at several temperatures from ~300 K to ~1000 K. 
The parameters obtained from these refinements will be 
described and discussed in the following section. 


IV. DISCUSSION 

The rhombohedral-to-cubic phase transition in GeTe 
is characterized by pronounced changes in the positional 
parameters rrce and :rTe- The irce (a-’Te) increases (de¬ 
creases) with increasing temperature. The initial increase 
(decrease) in XQ e (xx e ) is almost linear up to ~ 500 K and 
becomes non-linear while approaching the transition tem¬ 
perature T c = 600 K. One observes sudden raise (drop) 
of XQ e (xxe) at T c . Above T c , the XQ e and cctc take on 
the values characteristic for the cubic GeTe phase. These 
changes are also revealed by Aa; and Aa, see Fig. [4j as the 
























































































4 



300 400 500 600 700 800 900 


T(K) 

FIG. 4: (Color online) Temperature variations of (a) distor¬ 
tion parameter Ax and (b) deviation Aa of the distortion 
angle from 90°. The present NPD data (circles) are com¬ 
pared to the neutron diffraction experiments performed on 
single crystals by Chattopadhyay et al. [8] (squares). 

The continuous lines are just guides to the eye. 

distortion parameters are directly related to XQ e and XTe- 
Both Aa; and A a decrease continuously from Aa; = 0.14 
and Aa = 1.98° at 300 K to zero at T c = 600 K. The 
Ax and Aa obtained for single crystals of GeTe [8] show 
very similar behavior as that observed for powder sam¬ 
ples. They also progressively diminish with temperature, 
however with smaller slopes than those determined for 
our powder samples, and finally they approach zero val¬ 
ues at 700 K. The Aa; = 0 and Aa = 0 indicate that the 
R3m structure undergoes transformation into the Fm3m 
structure at T c . Here we note the difference in the transi¬ 
tion temperature T c between our powder samples and the 
single crystal GeTe samples investigated by Chattopad¬ 
hyay et al. |2]. Such a significant spread in transition 
temperature of GeTe compound has been known and is 
usually attributed to the sample stoichiometry as well as 
the free charge carrier concentration [251 . 

Distortions Ax and Aa has been considered respec¬ 
tively as the primary ( Q ) and secondary (e) order pa¬ 


rameters of the phase transition in GeTe [2] ■ The contin¬ 
uous shift of the Ge and Te atomic positions Ax which 
is directly related to the amplitude of the soft-phonon 
mode rj (Q) [IT] breaks the symmetry elements cor¬ 
responding to the four-fold rotation axis and a mirror 
plane perpendicular to this axis. It also destroys the 
symmetry inversion center of the high temperature cubic 
phase and induces polarization along the three-fold ro¬ 
tation axis. Thus, the primary order parameter Aa; ( Q ) 
is the driving force of ferroelectricity appearing in GeTe 
below T c . On the other hand, the Aa is actually the 
lattice strain e which breaks only the symmetry elements 
corresponding to the four-fold rotation axis and a mirror 
plane perpendicular to it, but it does not affect the sym¬ 
metry inversion center of the GeTe system. The neutron 
diffraction experiments on both powder and single crystal 
samples show the linear coupling between Aa and Aa; 2 , 
which conforms to the Landau theory [24]. The single 
crystal data [2] display, however, slightly smaller slope 
(0.89 x 10 -4 deg -1 ) than the present NPD data for pow¬ 
der samples (1.0 x 10 -4 deg -1 ). On the other hand, our 
theoretical calculations give the linear coupling between 
Aa and Aa; 2 of 1.2 x 10 -4 deg -1 . 

The short (s) and long (l) bond lengths determined 
from our neutron measurements on powder samples 
amount to s = 2.82 A and l = 3.18 A at room tem¬ 
perature. They are very close to those obtained from the 
neutron diffraction studies on single crystal samples [2]. 
The s and l bonds in the low temperature rhombohe- 
dral phase follow the course of XQe and XTe- They vary 
smoothly with temperature to reach the average value of 
~ 3 A above T c = 600 K, as shown in Fig. [E] This unique 
distance is certainly the Ge-Te bond length in the cubic 
phase which subsequently slightly grows at still higher 
temperatures due to the thermal expansion of the cubic 
GeTe lattice. Nevertheless the visible shift between the 
single crystal neutron diffraction data and those mea¬ 
sured on powder samples which arises from a difference 
in the respective transition temperatures, there is a good 
qualitative agreement between these two sets of data. 

Such a behavior as revealed by the Ge-Te bond lenghts 
in both single crystal and powder diffraction studies is 
however claimed m to be the case only in the so-called 
average structure. The PDF analysis of the total diffrac¬ 
tion data on GeTe shows that the high temperature phase 
exhibits two distinct bond lengths which hardly change 
across the phase transition m- The local scale structural 
distortions, evidenced by unequal Ge-Te bond distances 
above T c , suggested that the transition in GeTe could be 
of the order-disorder type nuns]. To obtain additional 
information about the local structure of GeTe above T c , 
we have simulated the PDF spectra of its high tempera¬ 
ture phase with the static structural lattice distortions 
generated according to the triply degenerate unstable 
soft-phonon mode of TJ symmetry m- This mode, 
while frozen, leads to relative displacements of the Ge 
and Te sublattices along the cubic cell diagonal. Results 
of these calculations for perfectly ordered ( Axq = 0) and 
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FIG. 5: (Color online) Temperature dependences of short (s) 
and long ( l ) Ge-Te bonds determined in the present NPD 
measurements (circles) and those provided by the single crys¬ 
tal neutron diffraction [SJ (squares). The continuous lines are 
just guides to the eye. 

disordered (Axi = 7.94 x 10~ 3 and Ax 2 = 14.56 x 10~ 3 ) 
high temperature cubic phases of GeTe are shown in 
Fig.© The distortions Aau and Ax 2 correspond re¬ 
spectively to the Ge-Te bond lengths of (si = 2.91 A, 
l\ = 3.10 A) and (s 2 = 2.84 A, I 2 = 3.19 A). Here we 
note that Ai 2 also reflects the Ge-Te bond lengths deter¬ 
mined above T c by using the experimental PDF data m- 
We observe that the simulated spectra of cubic GeTe with 
and without static distortions are essentially identical. 
Moreover, they are closely related to the experimental 
PDF data reported by Matsunaga et al. |T3] . This find¬ 
ing enables us to suggest that the PDF results can hardly 
provide a definite answer about the displacive or order- 
disorder type of the phase transition in GeTe compound 
since they probe the average static lattice distortions but 
not a dynamical nature of this transition connected with 
the phonon dynamics El- 

Figure [7] indicates that both powder and single crys¬ 
tal samples of GeTe exhibit temperature-induced vol¬ 
ume reduction AV at the ferroelectric rhombohedral-to- 
cubic phase transition. The unit cell volume decreases 
by AV ~ 0.6% at T c = 600 K, as indicated by our NPD 
data. The volume anomaly at the ferroelectric structural 
transition in GeTe can be considered as an electrostric- 
tion effect in analogy with magnetostriction effects at 
the magnetic ordering. Also, the volume contraction at 
T c has been interpreted by Chattopadhyay et al. [8j as 
the extra volume due to the presence of lone pairs in 
the rhombohedral phase and the absence of this excess 
volume in the high temperature cubic phase. Results 
of the DFT and QHA calculations performed for both 
the low temperature rhombohedral and high tempera¬ 
ture cubic structures of GeTe support to some extent the 
present experimental data. The agreement between cal¬ 
culations and experiment remains reasonable, except the 



5 10 15 20 25 30 

r (A) 


FIG. 6: (Color online) Pair distribution function G(r) cal¬ 
culated for the cubic GeTe phase with the static distortions 
Ax 0 = 0, An = 7.94 x 1(T 3 , and Ax 2 - 14.56 x 10~ 3 . 
Simulation are performed at the experimental r-range and 
Qma x = 15 A" 1 H3]. 
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FIG. 7: (Color online) Temperature dependence of the GeTe 
pseudo-cubic unit cell volume (F). The present neutron pow¬ 
der diffraction results (circles) are compared to the single crys¬ 
tal neutron diffraction data obtained by Chattopadhyay et al. 
[8]. Solid lines denote results of the DFT-QHA calculations 
for the low and high temperature phases of GeTe. 


close vicinity of T c , where the QHA approach is unable to 
reproduce the volume reduction at the onset of the phase 
transformation. Our calculations indicate that the vol¬ 
ume of the rhombohedral phase is higher than the volume 
of the cubic phase. Indeed, in the rhombohedral struc¬ 
ture the Ge and Te atoms are displaced from each other 
with respect to their position in the cubic structure. The 
relative shift of the Ge and Te sublattices along the body 
cell diagonal ( Ax ) is accompanied by angular distortion 
of the GeTe lattice (Aa) as well. Both distortions re¬ 
sult in the larger volume of the rhombohedral structure 
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in comparison with the rocksalt one. These distortions 
are dynamical in their origin as they are driven by the 
soft phonon mode , as shown and discussed in our 
recent paper m ■ The DFT calculations show that the 
rhombohedrally distorted GeTe is more energetically sta¬ 
ble than the undistorted cubic one. At the ground state 
the difference in their Helmholtz free energies amounts 
to 26 meV per formula unit, i.e., it lies in the range of 
thermal excitations. The present experiments also show 
that each of GeTe phases expands upon heating outside 
the temperature range where the phase transition occurs. 
The volume thermal expansion coefficient of the rhombo- 
hedral GeTe equals 4.59 x 10 -5 K -1 at 300 K, whereas 
it amounts to 7.67 x 10 -5 K _1 at 650 K for the cubic 
GeTe. 

The phase-change material GeTe is not a unique sys¬ 
tem undergoing temperature-induced volume collapse at 
the phase transition. The well-known example is the 
volume contraction of ice at its melting temperature. 
This phenomenon remains, however, unexplained quanti¬ 
tatively so far. The temperature-induced volume reduc¬ 
tion is also found to exist in solid-to-solid phase transi¬ 
tions in diverse condensed matter systems. Even earlier- 
known industrially important canonical ferroelectric ma¬ 
terial like BaTiC >3 shows such volume decrease at its 
ferroelectric phase transition at 393 K [23- One no¬ 
table example is the orbital order-disorder transition or 
orbital melting in the strongly correlated electron sys¬ 
tem LaMnOa [25i - the recognized parent compound of 
colossal magnetoresistive manganites that lure a host of 
condensed matter scientists for spintronics and other de¬ 
vice applications. Reduction of the LaMn 03 volume at 
the orbital order-disorder transition is assisted by the 
considerable change in the atomic mean-squared vibra¬ 
tional amplitudes [26] . Very similar effect is observed in 
the present studies on GeTe too. The isotropic temper¬ 
ature factors Bee and Bxe depicted in Fig. [8] increase 
with increasing temperature and show A-type behavior 
in the close vicinity of the rhombohedral-to-cubic phase 
transition. The Boe and Btc measured for single crys¬ 
tals of GeTe [8] are much more scattered compared to 
the present data. They are also limited to temperatures 
not exceeding T c and therefore they do not reveal such a 
characteristic change at T c as i?Ge and Btc in our exper¬ 
iments. 

It is well know that the atomic temperature factors 
remain sensitive to the immediate atomic environment. 
Experimental atomic thermal displacements contain both 
static and dynamic effects. The dynamic contribution 
to the atomic thermal displacements is gained from the 
present DFT calculations. The Uij tensor for both Ge 
and Te atoms in the R3m structure contains two indepen¬ 
dent components, namely U xx = U yy and U zz which rep¬ 
resent the atomic thermal vibrations perpendicular and 
parallel to the three-fold rotational axis, respectively (cf. 
Fig. 0- The site symmetries of the Ge and Te atoms in 
the FmSm structure constrain their mean-squared dis¬ 
placements to be isotropic and hence described by a sin- 



FIG. 8: (Color online) Experimental (symbols) and theoreti¬ 
cal (lines) isotropic temperature factors B determined for Ge 
(circles) and Te (squares) atoms in the rhombohedral and cu¬ 
bic phases of GeTe. The solid lines denote calculated B fac¬ 
tors for the rhombohedral and ideal cubic GeTe structures. 
Dashed lines correspond to calculated B factors for the cu¬ 
bic GeTe structure with static distortion Ax = 14.56 x 10 -3 
(s = 2.84 A, l = 3.19 A). Calculations are performed within 
the harmonic approximation. Theoretical isotropic tempera¬ 
ture factors are evaluated as B = 8n 2 ( U ), where {U) is the 
trace of Uij tensor [20] . 


gle parameter. Our calculations show that the thermal 
motions of Ge and Te atoms in the rhombohedral GeTe 
are almost isotropic as U xx ss U^ z and U xx « t/^ e . 
Similar dependences were also observed in the neutron 
diffraction study on single crystals of GeTe (SJ - The 
nearly isotropic Ujj tensor arises from very small differ¬ 
ence between the atomic on-site force constants along 
and perpendicular to the three-fold rotational axis in the 
rhombohedral GeTe. This negligible difference remains 
in accordance with the structural features of the rhom¬ 
bohedral phase for which the structural distortions of the 
Peierls-type are indeed small. The mean-squared vibra¬ 
tional amplitudes are obviously different for the Ge and 
Te atoms due to the difference in their masses. Also, they 
exhibit typical increase with increasing temperature ac¬ 
cording to the applied harmonic approximation with the 
Uij growing their values over two times between 300 and 
600 K. The harmonic approach is, however, unable to de¬ 
scribe properly the anomalous behavior of Uij in the close 
vicinity of T c . Nevertheless, it occurred to be helpful in 
explaining the behavior of atomic thermal vibrations in 
the cubic phase of GeTe. We notice that the slope of Bxe 
practically does not change it course with increasing tem¬ 
perature, except the tiny range of temperatures around 
T c = 600 K when it shows small discontinuity. There is, 
however, considerable jump in the thermal vibrations of 
the Ge atoms in the cubic GeTe phase as well as visi¬ 
ble change in the slope of Bq b above T c . The results of 
our calculations show the the on-site force constants of 
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Ge atoms in the cubic GeTe are almost two times lower 
compared to the respective on-site force constants in the 
rhombohedral GeTe structure, i.e., the Ge atoms in the 
cubic phase are more loosely bound inside the lattice 
which in turn allows for the larger values of the ampli¬ 
tudes of their mean-squared displacements. The mean¬ 
ingful decrease in the on-site force constants of the Ge 
atoms while going from the rhombohedral to cubic GeTe 
accounts for the increase in Bc e slope above T c = 600 K. 
It is interesting to note that the cubic phase of GeTe with 
incorporated static lattice distortions of Peierls-type Ax 
shows higher B-factors in comparison with the respec¬ 
tive factors for the ideal cubic structure (Ax = 0). In 
addition, the effect of Ax > 0 is much more pronounced 
for Bq b than for B^ e . Our calculations indicate that 
progressive increase of Ax results in the gradual growth 
and small changes in slopes of B-factors due to modi¬ 
fied force constants which decrease with decreasing Ax. 
In principle, this observation allows to distinguish the 
high temperature GeTe structure with persisting static 
lattice distortions (retained distinct short and long Ge¬ 
Te bonds) above the phase transition from that with the 
local distortions vanished (unique Ge-Te bond length). 

V. SUMMARY AND CONCLUSIONS 

Chattopadhyay et al. [5] carried out high tempera¬ 
ture single crystal neutron diffraction investigations on 
GeTe almost three decades ago when the powder neutron 
diffraction technique was not very mature. The present 
NPD investigations on GeTe have been carried out on a 
modern high resolution powder diffractometer at a spal¬ 
lation neutron source and have benefited from the enor¬ 
mous progress made during recent years in the neutron 
powder diffraction technique. The high temperature sin¬ 
gle crystal neutron diffraction is very time consuming and 
requires temperature stability for long periods whereas 
neutron powder diffraction is relatively free from these 
constraints. The extinction effects can be very large in 
single crystal diffraction whereas it is often practically 
negligible in powder diffraction. Also, the NPD technique 
appears much attractive in describing the phase transi¬ 
tions in relatively simple structures like GeTe due to its 
less sensitivity to the crystal domain structure which sig¬ 
nificantly complicates both the data collection and the 
data treatment. 

The results of our NPD measurements on GeTe are 
based on Rietveld refinement involving only Bragg in¬ 
tensities. An information contained in the background 
diffuse scattering from the sample is not taken into ac¬ 
count, and hence our results are related to the average 
structure, but not to the local structure or dynamics. 
Although the current experimental results could not un¬ 
ambiguously resolve controversy about the nature of the 
phase transition in GeTe, viz. whether this phase transi¬ 
tion is of displace or order-disorder type, they allowed for 


more accurate determination of the temperature depen¬ 
dences of the lattice and structural parameters of GeTe 
compared to the previous studies [5|. The measured vari¬ 
ation of the GeTe volume over a wide temperature range 
is however a robust result of Bragg diffraction and the 
volume anomalous behavior at the phase transition is the 
most important result. The temperature variation of the 
structural parameters across the rhombohedral-to-cubic 
phase transition provides us microscopic mechanism be¬ 
hind the volume discontinuity and enables us to remove 
the veil of mystery around this transition. 

It is interesting to note that the interpretation of the 
local probes results urns] viz. the PDF analysis of the 
total diffraction intensities or EXAFS, give an impres¬ 
sion that apart from a small linear thermal expansion 
nothing happens during the phase transition, namely the 
short and long Ge-Te bond lengths do not change within 
experimental errors in the whole wide temperature range 
investigated (300-800 K). One also wonders whether a 
hypothetical static model compatible with the bond dis¬ 
tances obtained by the local probes would reproduce the 
robust result of volume contraction at the rhombohedral- 
to-cubic phase transition in GeTe. Perhaps some kind 
of a dynamical model, though definitely difficult to con¬ 
struct, might reproduce the volume decrease. Hence, an 
answer to the question whether transition is of displacive 
or order-disorder type lies in the dynamics of phase tran¬ 
sition and the role of soft modes, as has been recently 
shown by the DFT and phonon calculations m The 
conventional PDF analysis of the total scattering data 
does not analyze the energy and therefore it lacks an in¬ 
formation about dynamics of the phase transition. The 
high temperature phase contains the dynamics of broken 
symmetry phase and the PDF local probe just sees the 
snap-shot of low frequency soft phonon mode that hap¬ 
pened to enter into the window of local probe. The two 
distinct Ge-Te bond lengths seen by the local probes at 
the high temperature phase are just dynamical signature 
of the low temperature phase still persisting at temper¬ 
atures exceeding T c . Since the PDF measures instan¬ 
taneous structure and cannot distinguish between static 
and dynamic correlations, interpretation of the PDF re¬ 
sults based solely on the static bond distances could be 
verified. 
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